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Abstract 

We propose a new model of cluster growth according to which the probability that a 
new unit is placed in a point at a distance r from the city center is a Gaussian with 
mean equal to the cluster radius and variance proportional to the mean, modulated 
by the local density p{r). The model is analytically solvable in d = 2 dimensions, 
where the density profile varies as a complementary error function. The model re- 
produces experimental observations relative to the morphology of cities, determined 
via an original analysis of digital maps with a very high spatial resolution, and helps 
understanding the emergence of vehicular traffic. 
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1 Introduction 

Important problems related to cluster growth processes occur in a number 
of different disciplines, ranging from physics to biology and transportation 
engineering, and several microscopic models have been proposed to describe 
the growth of both compact clusters, like crystals or tumors [2,1], and fractal 
clusters, like colloidal aggregates or snowflakes [3]. All of these models are 
characterized by the presence of an 'active' zone on the surface of the cluster 
where the cluster growth takes place. Despite the simplicity which charac- 
terizes the microscopic dynamics of these cluster growth models, analytical 
solutions for the temporal evolution and for the spatial dependency of cluster 
properties, like the density profile or the width of the active zone, are difficult 
to obtain. To this end one usually resorts to models for the evolution of the 
cluster surface, such as the Kardar-Parisi-Zhang model [4], or to extensive 
numerical simulations [5,6,7]. 
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These numerical simulations have suggested that, at least in the case of the 
Eden Model and of the DLA model, the radially averaged probability P(r, N)dr 
that the (A^ + l)th cluster unit is deposited within a shell of width dr at a 
distance r from the center of mass of the cluster is well approximated for r 
and N large by a Gaussian distribution. 
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with mean r^v oc and variance a oc N'^'; the scaling exponents u and u' 
are model dependent. This growth probability distribution leads to a density 
profile p(r, N) of a cluster of size N in d dimensions, 

1 ^ 

pir,N)^^-^Jp{r,N')dN', (2) 



where Sa is the surface area of the d dimensional unit sphere {S2 = 2ti, S3 = tt), 
which can be evaluated in the limit N ^ 00 and r large [5]: p(r, N — 00) oc 

j.—d+l/v 

Starting from Eq. 1 and from some considerations about the asymptotic be- 
havior of P{r, N) in the r — > and r — > 00 limits, in this paper we elaborate a 
new model for the growth of compact clusters, the random walk growth model 
(RWG). This is based on the simple idea that the cluster radius grows as a 
random walker subject to a drift, which gives a growing probability distribu- 
tion Prwg{^-,N) oc r'^~^P{r,N). We solve the model in d = 2 dimensions, 
showing that the density profile varies as a complementary error function. As 
an application of the proposed model we have studied the morphology of sev- 
eral European cities, which are growing clusters. Via an original analysis of 
digital maps with a very high resolution [8] we have determined the spatial 
dependence of their density of streets Ps(r), which appears to be very well 
described by the RWG model. 



2 The Random Walk Cluster Growth Model 



In a large number of cluster growth models (DLA, Eden, Sohd-on-sohd, Ran- 
dom Deposition, . . .), a cluster grows as a new cluster unit is placed near an 
existing one. Therefore, in order for a cluster to grow in a given location r, 
at least a cluster unit must be present near f. In this respect, it is suprising 
that the growth probability of Eq. 1 depends on the cluster density p(r, A^) 
only through r^v and cr^v; instead, one would have expected the radially aver- 
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aged probability of placing a cluster unit in a shell at a distance r from the 
cluster center to be proportional to the number of cluster units which occupy 
the shell, i.e. P{r,N) oc p{r,Ny~^ oc r^f~^, where Df = l/u is the fractal 
dimension of the cluster. 

We want also to point out that if P{r,N) oc r^f~^^ then Eq. 2 predicts the 
cluster density to diverge as p(r) oc r^f~'^ when r — 0. We therefore expect 
a crossover in P(r, N) which must be proportional to r'^~^ when r — > 0, and 
proportional to r^f~^ when r ^ oo. 

Inspired by Eq. 1 and keeping in mind the above considerations, here we 
define a new model for the growth of compact clusters, where no crossover 
in the radial growing probability is expected a.s Df = d, which gives rise to 
a growing probability distribution PRWG{f-,N) oc r'^~^P{r,N). The model is 
defined by assuming 1) that the cluster mass N is related to the mean cluster 
radius rjv by = {tm/tqY in d spatial dimensions; and 2) that the radius 
of the cluster evolves as a random walker subject to a drift: at each updating 
step the radius varies of a quantity taken from a distribution with mean w > 
(drift velocity) and variance u^. 

Under these assumptions the radially averaged probability -Prvkg(^) N)dr that 
the (iV+ l)th cluster unit is deposited within a shell of width dr at a distance 
r from the center of mass of the cluster is 

^RWG(^^) = ^^(^^)' (3) 



where P(r, A^) is given in Eq. 1, //^ = drr^P{r, N) and we assume k — d—1 

as discussed above. The growth dynamics sets u = 1/d and v' — v/2 = l/2d. 
Note that for A; < — 1 the cluster density diverges when r — 0, while for 
k > d — 1 it does not decrease monotonically, and has a maximum at r > 0. 

We restrict our analysis to d = 2 dimensions, which is the more interesting 
case from a physical viewpoint, but an analytical treatment is also possible 
for d — 3. 

According to Eqs. 2 and 3 the density of a cluster of size N — {rf^/roY, is 
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r (arbitrary units) 

Fig. 1. Model of cluster growth. The probability -P|^"wg(''' -^)/^'^^^ '^^ building a 
cluster unit in a point at a distance r from the center is Gaussian centered around 
the cluster radius r^v = tqN^^'^ with variance af^ = a'^r^N^I'^ /v, Eq. 3. Here we plot 
in arbitrary units the cluster density p{r, N) (Eq. 3) and the corresponding growing 
probability -Pr;wg(^' fN = ^ (dashed lines) and = 24 (full lines) in 
d = 2 spatial dimensions. 

and via the change of variables y = n~^^'^"'^, 
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The last term of Eq. 5 can be safely neglected in the case of large clusters, 
^^Af/o'AT ^ 1 ^ , leading to the following approximation 



p(r, A^) ~ Poo + 



Pmax ~ Poo J. ( r 



-Erfc 



In 



Vy2a^ 



(6) 



according to which the cluster density decays as a complementary error func- 
tion. Here poo = p('", 0) is the density before the growth of the cluster, or 
equivalently the density far away from the cluster center (r ^ tat), and 
Pmax — Poo = (27rrof The maximum value of the density is achieved in the 
cluster center; for large clusters r^/aN = ^/r^v /a ^ 1 and p(0, iV) ~ pmax- 
A density profile similar to that of Eq. 6 is given by the solid-on-solid growth 
model in the high deposition limit [2]. 

^ A numerical evaluation of the usual L2 norm of p{r, t) shows that an overestimate 
of the relative error made neglecting the last term in the case of real clusters (see 
Table 1) is 10"^, which is obtained with A = cr^/rjv = 0.6 Km and r^r = 1 Km. 
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Fig. 2. The density of streets ps{x,y) of the city of Paris (upper panel) and of the 
city of London (lower panel). Each cell of this picture covers an area of size 1.9 x 1.7 
Km^ and indicates, according to the color code showed on the right, the percentage 
of the cell surface which corresponds to a street. The center of the city {x = y = 0) 
is located on the lie de la cite for the case of Paris, and near Westminster for the 
case of London. 

A schematic drawing explaining the growing dynamics of a cluster according 
to the proposed model is shown in Fig. 1. As a cluster grows A^, tn and <jn 
increase; -Pr;wq {t, N) shifts versus higher values of r and becomes shorter 
and wider. 

We do expect this model to be able to describe the growth of compact clusters, 
characterized by fluctuating and growing interface. This is the case of real 
cities, which grow as new buildings are constructed, usually in the suburbs, 
as we show in the following section. 



3 Empirical validation of the model 

We have validated the model by studying the radial dependence of the density 
of streets Ps(r) of several European cities (here the subindex s is used to 
indicate quantities related to streets), which has been obtained by analyzing 
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City 


TN (K-m) 


A (Kmj 




Poo 


Paris 


Li .( 


U.oi 


U.OD 


U.Ud 


London 


14.5 


0.54 


0.36 


0.05 


Rome 


8.7 


0.31 


0.20 


0.07 


Modena 


1.2 


0.34 


0.16 


0.05 



Table 1 



Parameters characterizing the radial dependence of the density of streets of various 
European cities, obtained by making a fit of empirical data to Eq. 5. All reported 
values have relative errors smaller than 5%. 

digital maps [8] with a very high spatial resolution (1 pixel correspond to 
around 3.3 m^). The density of street of a given area A is given by — 

jn^^ where is the total number of pixels which cover the are A, and 
is the number of pixels which cover the area A, and also lie on a street. As far 
as we know this is the first empirical determination of the radial distribution 
of the density of streets, and no previous studies achieve such a high resolution 
in the determination of the morphology of a city. 

As an example of the result of our digital image analysis we show in Fig. 2 
the density of streets Ps{x,y) of the cities of Paris, Prance (upper panel) and 
London. U.K. (lower panel). Por visualization purposes the images have a 
resolution 1.9 x 1.7 Km^, smaller than the actual resolution reached in our 
study. As expected the density of streets reaches higher values in the city 
center (x ~ j/ ^ 0), and smaller values in the suburbs. This is quantified in 
Pig. 3 where we show the radial dependence of the density of streets Ps{r) 
of several European cities (London (U.K.), Paris (Prance), Rome (Italy) and 
Modena (Italy)) which can be considered to a first approximation round. Ps{r) 
appears to be always well fitted by the functional form of Eq. 5. 

Prom the fitting parameters one can compute the radius of the city rjv, the 
variance a]^ of the growing distribution or a characteristic growing length A = 
(y%/r^ = (j'^/v, and the values of the density of streets in the city center /Omax 
and in the suburbs Poo- These values are reported in Table 1. The parameter 
rjv is a measure of the radius of a city, and therefore varies from city to city. 
The same is true for Pmax^ Paris and London, which are characterized by the 
presence of large streets within the city center, have a higher value of pmax than 
Rome and Modena. The parameter po is the same for all considered cities, and 
therefore appears as a unifying property of the street network of the suburbs 
of all cities. Pinally, the parameter A assumes similar values for the cities of 
Paris, Roma and Modena, and a higher value for the city of London. As A 
controls the variation of the variance of the distribution via a% — Ar^v we 
conclude that London is less compact than other European cities. 

The possibility of describing the radial dependence of the density of streets 
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Fig. 3. The density of streets as a function of the distance from the city center 
Ps{t) for the cities of Modena, Rome, Paris and London. This last curve is shifted 
verticaUy of 0.1 for clarity. Plain lines are fit to Eq. 5. 

of all considered cities with the same functional form (Eq. 5) suggests the 
existence of a common mechanism underlying the growth of all cities, which 
appears to be captured by the RWG cluster growth model. The success of 
Eq. 5 in explaining the empirical data questions the common assumption of 
exponential decay of the population density [10,12,11]. In fact, (at least for 
r < R) it is reasonable to assume a proportionality between density of streets 
and population density (or density of urbanized area). 

Our model reproduces the radially averaged properties of a city. In principle 
the radial dependency of Eq. 5 could be used in order to reproduce spatial 
correlations in urban settlements, as in [12]. 



4 Communication and congestion 

In many physical situations a cluster plays the role of a transportation infras- 
tructure, as it is the backbone over which something moves. This is the case, 
for instance, of communication or transportation networks, both biological 
and artificial, which clusterize near 'vital' centers. In these cases the topology 
of the cluster influences the dynamics of the transportation processes which 
occur over it [9]. Particularly, large spatial gradients in the cluster area may 
lead to congestion and to a dramatic slowdown of the dynamics. Here we dis- 
cuss these influences in the case of clusters with a density given by the RWG 
model (Eq. 5). We consider for definiteness the case of vehicular traffic over 
the street network of a city, but the results are more general. 

Empirical evidences shows that there is a transition from free traffic flow to 
congested traffic flow as the car density pcar increases. For instance highway 
traffic becomes congested when pcar > Pcar^ — 30 vehicles/Km [13]. The knowl- 
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Fig. 4. Area occupied by the streets in a circular annulus of width A = 1 Km 
centered on the city center, S{r) = 2irrAps{r), as a function of the distance from 
the city center r. The curve refers to the city of London, but similar curves are 
found for the cities of Paris, Rome and Modena. The non-monotonic behavior of 
S{r) influences vehicular traffic: cars moving outward are compressed in the region 
17 < r < 36 Km where S{r) decreases. 

edge of the density of streets allows for an estimate of the car density. Given 
cars in a region A one can estimate pcar = PN/Sa where Sa = Ja Ps{x, y)dxdy 
is the area occupied by the streets in the region A, and P a coefficient of 
proportionality. 

Let's consider what happens when cars, initially located in a circular an- 
nulus of radius r and width A around the city, move outward. As their mean 
distance from the city center r increases the car density varies as Pcari^) = 
f3N/S{r), where 



r+A 



S{r) = 2-K 'rjps{r')dr'2TcrA 



Poo + 



Pmax — Poo J. ( r 



-Erfc 



(7) 



is the area occupied by the streets in the circular annulus we are considering. 
The car density is higher where S{r) is smaller, and increases in the regions in 
which dS/dr < 0. These regions are those which most probably act as traffic 
bottlenecks. 

For the city of London, Figure 4 shows that S{r) increases for r < 17 Km 
and for r > 36 Km, when the density of streets varies very slowly. In the 
range 17 < r < 36 Km of fast variation of the density of streets S{r) de- 
creases. It follows that in this range of distances from the city center, i.e. at 
the city boundary, traffic is maximum. This is therefore the area where new 
transportation infrastructures and strategies for traffic optimization should 
be concentrated. The non-monotonic behavior of S{r) also characterize Paris, 
Rome and Modena. 
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5 Conclusions 



In this paper we have introduced and studied a new model of cluster growth. 
This model is inspired by numerical results relative to the DLA and to the 
Eden model, but it takes into account the fact that the probability of building 
a cluster unit in a given point depends on the local density. In this model 
the probability of building a cluster unit in a point at a distance r from 
the cluster center is a Gaussian centered around the radius of the cluster, 
with variance proportional to the radius, modulated by the local density p{r). 
An explicit solution of the cluster properties in d = 2 dimensions, which is 
not available for other cluster growth models, shows that the cluster density 
decays as complementary error function. We have validated the model via a 
high resolution study of the density of streets of several European cities, and 
we have discussed the relation between the cluster topology and the dynamics 
over the cluster. 
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